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In this Brief Report, a non-standard solution to the solar neutrino problem is revisited. This so- 
lution assumes that neutrino flavors could have different couplings to gravity, hence, the equivalence 
principle is violated in this mechanism. The gravity induced mixing has the potential of accounting 
for the current solar neutrino data from several experiments even for massless neutrinos. We fit 
this solution to the total rate of neutrino events in the SuperKamiokande detector together with 
the total rate from other detectors and also with the most recent results of the SuperKamiokande 
results for the recoil-electron spectrum. 



I. INTRODUCTION 



At present several neutrino detectors indicate that there is a major discrepancy between the Standard Solar 
Model (SSM) predictions for the solar neutrino flux 111 and its observed value. Recently, five neutrino experiments 
(GALLEX SAGE j|, Homestake ||, Kamiokande [3 and SuperKamiokande |^|) reported that the observed neu- 
trino flux is less than the values predicted by the SSM for different neutrino energies. The first three experiments 
detect the solar electron neutrinos, u e , via absorption processes so they do not give any information about the direction 
and energy spectrum of the observed neutrino events. The last two water Cherenkov detectors employ v e -e scattering 
so the direction and energy distribution of the detected neutrino can be fairly determined. This would allow us to test 
theoretical models not only against the measured total rate but also against the energy spectrum of the neutrino flux. 
The analysis of the first 504 days data from SuperKamiokande |t]](SK), when combined with the data from earlier 
experiments, provide us with important constraints on the famous MSW effect S and vacuum oscillation solutions of 
the solar neutrino problem. The MSW effect involves the standard neutrino interactions with the matter background 
in the Sun, leading to an enhanced v e -v x resonant transition. The "just so" vacuum oscillation solution explains the 
suppression using only vacuum oscillations of the neutrinos on their way to the Earth. This solution is in fact in 
slightly better agreement with the recent data from SK than the matter-induced MSW effect However, it requires 
some "fine tuning" of the Sun-Earth distance to account for the apparent deficit. Both of these solutions require that 
neutrinos have non-degenerate masses; they will be referred to as the Mass-Induced (MI) scenario. 

In this report, we explore yet another possible solution of the solar neutrino problem. However, it requires a rather 
unorthodox assumption that the neutrino flavors couple differently to gravity. This is a statement of the violation 
of the weak equivalence principle which stipulates that all matter couple equally to gravity. This solution has been 
suggested by different authors p0j , pi| ; here, we consider the Gravitational-Induced (GI) scenario in the light of the 
most recent results of the total rate and recoil-electron spectral shape from the SK detector. It should be noted that 
the pure GI mechanism does not require the introduction of neutrino masses and thus can be consistent with a zero 
or degenerate mass scheme for light neutrinos (v e , or v^). We use a two-flavor analysis treatment in our analysis for 
simplicity. A complete GI three-flavor analysis treatment was studied ]T^] ; however, the present experimental results 
are not enough to constrain the parameter space for a three-flavor analysis. So the assumption made here is that the 
GI mechanism, if it were to be the cause of the solar neutrino problem, is due to only two neutrino flavors. A mixed 
GI+MI scenario where neutrinos are assumed to be massive was studied by the authors of |l3| assuming that the MI 
effect is the main cause of the neutrino flux deficit. We do not cover this possibility here. 

In addition to the GI scenario considered here, there are other scenarios which are related to the violation of 
equivalence principle. The first involves the neutrino coupling to the massless dilaton field which arises from String 
Theory pUT^ ]. This scenario allows for a vacuum oscillation solution of the solar neutrino problem, but there is no 
clear distinction to be made in the present neutrino experiments between it and the conventional mass mixing. The 
second one proposed by Glashow et al. |]l6| assumes that Lorentz invariance is violated during neutrino propagation. 
The analysis is equivalent to the GI scenario with a redefinition of parameters. In this scenario, different neutrino 
flavors are allowed to have different velocities under the assumption that the gravitational potential does not 
change appreciably over the distance of interest. 

The article is organized as follows. Section II reviews the 2-flavor formalism for the pure GI mechanism of neutrino 
oscillations, detailing the survival probability of the neutrinos during their propagation in a gravitational potential. 



1 



Section III presents a x 2 analysis of the solar neutrino data which is used to constrain the relevant GI parameters. 
Some concluding remarks are given in Section IV. 



II. GI FORMALISM FOR TWO NEUTRINO FLAVORS 



The violation of the equivalence principle has been tested in several famous experiments by Eotvos, Dicke et al. , 
etc. |T?j] . The current bounds on the extent of the violation still allow an appreciable parameter space for the GI 
mechanism to have an effect on the neutrino propagation in regions where we have large gravitational fields. We 
discuss the case of a pure GI mechanism where the neutrinos are assumed to be massless. The neutrino wavefunction 
tjj l (x), where i denotes its flavor, obeys the gravity modified Dirac equation, 

[r,^ + h^(x)]d fl d u ^ i (x) = 0, (2.1) 

where (x) are the perturbations in the metric due to the gravitational potential <fr(x). The h^'s can be 
parametrized as fiG^, where fi are the flavor coupling constants to the gravitational field G^ v ' . The constants 
fi can be regarded as post-Newtonian (PPN) parameters. If Einstein's equivalence principle is not obeyed, the 
parameters fi will not be equal for different neutrino flavors. 

In the GI mechanism the neutrino weak basis v a does not coincide with the gravity basis v a and the neutrino 
eigenstates mix according to 

"a - U^ a , (2.2) 

where U G is a unitary matrix depending on the gravity mixing angle 9q. The correspondence between the GI and 
the MI formalisms is obvious if we perform the following replacements: 

@G — ► &MI 



and 

2E6f$^— . (2.3) 

Sf = f e — f x denotes the difference between the gravitational couplings of v e and v x where we will be considering 
the transition between the electron-neutrino flavor v e and another flavor v x (z^, or v r ) which is relevant to the solar 
neutrino problem. Sm 2 is the corresponding squared mass difference in the MI effect. 

In the presence of ordinary matter, the propagation equation for massless neutrinos in a gravitational field can be 
written as 

££)(::)■ <"> 

where E is the neutrino energy and 6j = \&\ 5f. We shall use £7 in our analysis instead of Sf due to the ambiguity 
in choosing the potential $(r) jnjfj]. ASj = \/2GfN £ /2E represents the charged weak interaction term between the 
neutrinos and the electrons in matter with a number density N e . In vacuum, the survival probability can be easily 
obtained, 

ttI 

P(v e ^v e ) = l-sin 2 28 G sm 2 —, (2.5) 



where Ac is the GI oscillation length defined as 



XG = imtf- (2 - 6) 



1 Many authors suggest that the correct choice of the potential lies in choosing the potential of the Great Attractor, <£> ~ 10 -5 , 
not the solar potential, $ ~ 1CF 6 . We will assume the former potential to be the dominant in our analysis especially that it 
has no considerable variation over the relevant distance scale. 
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The energy dependence of X G is obviously different from the case of the MI oscillation length which is inversely 
proportional to E. This has led several authors (ll]] to suggest that the best way to distinguish between the two 
mechanisms is by measuring the neutrino energy spectrum. 

Now we turn to the generalization of Eq. (2.5) in the presence of matter. By looking at the propagation equation, 
Eq. (2.4), we note that we can have a resonant transition between v e and v x only if 5f is positive and 

2£J7 _^62. (2 . 7) 



cos 29 G 

The adiabaticity condition in the GI mechanism reads 

8-fE sin 2 29 



1 



1 dN e _ 1 ri$ 
N, dr $ dr 



cos 29 G 



» 1, (2.8) 



resonance 



where we assume in our calculation that the variation of the gravitational potential of the Great Attractor with 
distance is slow, i. e. A 4^ <C T^-^r^- Using the analogy to the derivation of the MSW survival probability, we can 
write the survival probability in the presence of matter as 

P(v e -> v e ) = i + Q - P C J cos 29 G cos 29 Gm , (2.9) 

where the matter mixing angle 9a m is defined as, tan29 Gm = sin 29 c4 / {cos 29 g — A) and P c is the level-crossing 
probability given by its standard form for an exponentially varying density, 

P c = (exp[-^7(l - tan 2 9 G )] - exp[-27T7 cos 29 G j sin 2 29 G \)j 

( 1 - exp [- 2tt7 cos 29 g / sin 2 29 G ] ) . 

(2.10) 

The survival probability is plotted versus Edj for different values of sin 2 29 G in Fig. [j]. 



III. NUMERICAL ANALYSIS 

In the following discussion, we present a x 2_ana lysis of the most recent solar neutrino data in the light of the 
GI mechanism. In Section Ilia, a x 2 analysis on the total rate of the neutrino flux is first performed to obtain 
the constraints on the GI parameters (Sj, sin 2 20(3) • The following section, Section Illb, deals with the analysis of 
the recoil-electron spectrum from the SK detector. Both analyses are then combined to constrain the relevant GI 
parameters. 



A. Total Rates 

The experimental results of the total rates for the four neutrino experiments used in this article are shown in Table |. 
We use the predictions of the Bahcall-Pinsonneault standard solar model of Ref. [jl8| (BP98) assuming the Institute 
of Nuclear Theory (INT) estimate []l9| for the B 8 production cross section. Our procedure takes into account the 
theoretical uncertainties in the theoretical estimates of the neutrino flux. The \ 2 function is given by J^] 

X? ates = E - R? v ) V^ 1 {Rf - Rf p ) , (3.1) 

ij=l,...,4 

where i?* h ^ cxp ' represents the GI theoretical predictions (experimental values) of the neutrino detection rate divided 
by the SSM predictions for the «-th experiment. Vij is the error matrix which contains the theoretical uncertainties 
as well as the experimental errors for each experiment. Fig. || shows the 90% and 95% C. L. regions for the GI 
parameters and sin 2 29 G with Xmin = 0-46. This is to be compared to our estimate of the MI best fit corresponding 
to the MSW small angle solution to the solar neutrino problem, where xLm = 0-62. In Fig. ||, we see that there are 
two allowed regions: the first is at 5^/ ~ 10 -18 and sin 2 29 G ~ 10 -3 (small angle solution); while the second region 
appears at a smaller value of £7 ~ 10 -21 and maximal mixing where 9 G < 7r/4 (large angle solution). The analogy 
between the standard MI solutions and the GI solutions is understood since the GI calculations involved with the rate 
analysis come straight from the MI calculations with the parameter replacement of Eq. (2.3). 
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B. Recoil electron spectrum 



We next consider the analysis of the SK recoil electron spectrum using the following % 2 : 

Xs 2 pect r a= E Hf-C P )^r 1 M h -C P )- ( 3 - 2 ) 
z,j=l,...,16 

Here 0* h ( exp ) is the i-th energy bin theoretical (experimental) value for the measured flux divided by the SSM 
prediction. The error matrix used here is given as V/j = af tat <jj tat 5ij + (j| ys cr^ ys ||. a is the flux normalization 
parameter which is allowed to vary independently of £7 and sin 2 20 g- This variation allow for the testing of measured 
spectrum and not the overall rate of SK. Fig. || shows the ratio of the GI predicted flux to the SSM predictions 
for the 8 B neutrino flux plotted against the recoil electron energy. The predicted ratio and the experimental results 
are normalized at 11.7 MeV. Fig. |^ presents the spectral shape analysis exclusion region at the x 2 = 25.0 level for 
14 degrees of freedom (16 energy bins minus two free parameters) together with the allowed regions from the rate 
analysis. We note that both the small angle solution and the large angle solution from the rate analysis are not 
affected. 

In the spectral shape analysis, we allow the hep flux contribution to vary as suggested in Ref. (2(]] to see whether 
we will get any effect on the exclusion region or not. It is found that an hep flux enhancement of about 10 3 is needed 
to have any appreciable effect on the exclusion region. In any case, the exclusion region does not intersect with any 
of our allowed regions from the rate analysis. Thus, the spectral shape analysis based on the recent SK data yields 
no additional constraints on the GI allowed parameter space. 

IV. DISCUSSION 

We have shown in this report that the GI scenario is still a viable solution to the solar neutrino problem. From just 
the rate analysis, we obtain two allowed regions in the GI parameter space: the first is at £7 ~ 10 -21 and maximal 
mixing, while the second is at a smaller value of 5j(~ 10~ 18 ) and sin 2 20 ~ 10 -3 . The spectral shape analysis has no 
effect on these allowed regions. We should also note that the combined best fit regions are still allowed by the current 
experimental bounds on the violation of the equivalence principle and by the neutrino accelerator experiments 
Although the allowed parameter regions for the GI scenario are quite small and may seem unlikely, they are by no 
means excluded by the available data. Thus, to obtain the final word about either the MI or GI explanations to the 
solar neutrino problem, more data are needed especially on the solar neutrino spectral shape from the SK and SNO 
detectors so that we can eventually accept or reject one of these scenarios. 
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TABLE I. The current experimental results and SSM predictions for the total rate in four neutrino detection experiments. 
The rate is measured in SNU for all of the experiments except the SuperKamiokande detector. The SSM predicted rates are 
based on the Bahcall-Pinsonneault standard solar model (BP98) of Ref. [18] . The theoretical errors are at the la level. 



Experiment 



Experimental Rate 



SSM predicted Rate 



Threshold Energy 



Homestake (CI) [§ 
GALLEX (Ga) @ 
SAGE (Ga) § 



2.56 ±0.16 ±0.14 




0.81 MeV 



77.5 ±6.2 



±4.3 
-4.7 




0.24 MeV 



66.6 



±7.8 
-8.1 




0.24 MeV 




2.44 ± ±0.05 



±0.09 
-0.07 



6.5 MeV 



FIG. 1. The survival probability P(v e — > v e ) in the GI scenario. It is plotted against the product of the neutrino energy E 
and the measure of the equivalence principle violation £7 (a = E8-y in units of 10 -12 eV). 

FIG. 2. Contour plot of the allowed regions in the (Sj, sin 2 26) plane. These are found from the results of the total rate 
for the four neutrino experiments (HomeStake, GALLEX, SAGE and SuperKamiokande). The GALLEX and SAGE results 
have been combined in the analysis. The solid curve represents the 95% C. L. region while the dashed one indicates the 90% 
C. L. region. The star shows the position of the best fit solution at Xmin- 

FIG. 3. The ratio of the GI predicted recoil-electron spectrum F(E e ) to the SSM spectrum Fssm(E £ ) for different values 
of 57 and 80 ■ The ratio is normalized at E e — 11.7 MeV. The error bars indicate the statistical and systematic errors of the 
experimental data. 

FIG. 4. Contour plot of the exclusion region in the (67, sin 2 28) plane, the dotted curve shown is an isocontour at \ 2 — 25.0 
for 14 D. O. F. This has been based on the recent results of the recoil-electron spectral shape from the 504 days SK data run. 
The regions allowed from the rate analysis are also shown. As in Fig. 1, the dashed curve represents the 90% C. L. while the 
solid curve represents the 95% C. L. 



2 The SuperKamiokande rate is measured in 10 6 cm 2 s 
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